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Abstract 

It is shown that the system of two coupled harmonic oscillators possesses many interesting 
symmetries. It is noted that the symmetry of a single oscillator is that of the three-parameter 
group Sp( 2). Thus two uncoupled oscillator exhibits a direct product of two Sp( 2) groups, 
with six parameters. The coupling can be achieved through a rotation in the two-dimensional 
space of two oscillator coordinates. The closure of the commutation relations for the gen- 
erators leads to the ten-parameter group Sp(4) which is locally isomorphic to the deSitter 
group 0(3,2). 


1 Introduction 

Since the classical mechanics of two coupled harmonic oscillators is known to every physicist, there 
is a tendency to believe that this oscillator problem is completely understood and that nothing 
new can be learned from it. We are writing this note because there are so many new lessons 
to learn from the system of coupled oscillators. The system shares symmetries with a number 
of physical models of current interest, such as the Lee model in quantum field theory [1], the 
Bogoliubov transformation in superconductivity [2, 3], two-mode squeezed states of light [4, 5, 6], 
the covariant harmonic oscillator model for the parton picture [7], and models in molecular physics 
[8]. There are also models of current interest in which one of the variables is not observed, including 
thermo-field dynamics [9], two-mode squeezed states [10, 11], the hadronic temperature [12]. 

From the mathematical point of view, the standard approach is to construct a suitable repre- 
sentation of the symmetry group after writing down its generators. The symmetry group in the 
present case is 5p(4) with ten generators [4, 6]. However, it is extremely difficult to study physics 
in terms of ten parameters. We should somehow start with a smaller number of parameters. 

For example, let us consider the three-dimensional rotation group with three generators. We 
need only two generators to describe rotations [13]. The third generator is produced during 
the process of constructing a closed set of commutation relations. For the coupled oscillators, a 
reasonable approach is to start with simpler groups describing two uncoupled oscillators. We can 
then introduce an additional generator to couple the two oscillators. The number of generators 
of the resulting group may be larger than the sum of those for the two starting groups plus the 
additional generator to couple them. 
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The process of constructing a larger group from two smaller groups is quite common in physics. 
We are quite familiar with the ’’direct product” and ”semi-direct product,” where the number of 
generators is the sum of those for the two smaller groups. We shall use the word “construction 
of the group by soldering two subgroups,” when the resulting group has more generators than all 
those for the starting groups and the additional generator(s) to couple them. We need this new 
word ’’soldering” in order to reduce the number of input parameters in physics. 

The soldering process takes different forms. We can construct the three dimensional rotation 
group by soldering the one-parameter rotation group around the x axis and another rotation group 
around the y axis. The resulting group has three generators. As we shall see in this paper, we 
can construct the 0(3,2) deSitter group by soldering two 0(3, 1) Lorentz groups. In this case, we 
solder them by observing that the two 0(3, 1) groups share the same rotation group. We start 
with nine generators, but the resulting 0(3,2) deSitter group has ten generators. 

Since the symmetry group of each uncoupled oscillator is the three-parameter Sp( 2) group, and 
since it is likely that one more group operation is needed to couple the system, w r e start here with 
seven generators. We shall see in this paper that the resulting group is Sp(4) with ten generators. 
It is easier to study physics with seven generators than with ten. 

It is also shown in this paper that the Sp( 4) symmetry does not exhaust all possible symmetries 
of the coupled oscillator system. It is noted that the group Sp( 4) is a subgroup of a lager group 
5Z/(4, r). Possible physical implications of this larger symmetry group are discussed. 

In Sec. 2, we shall study linear canonical transformations in the four-dimensional phase space 
consisting of two pairs of canonical variables. It is noted that the symmetry group is Sp( 4) which 
is locally isomorphic to the 0(3,2) deSitter group. In Sec. 3, we shall see how the 0(3,2) group 
can be constructed from two (3 + l)-dimensional Lorentz groups. In Sec. 3, we shall construct 
the symmetry group of two coupled oscillators from the symmetry group of each oscillator. 

Section 5 contains a new parametrization of the coupled oscillator system which is consistent 
with that of the symmetry group. In Sec. 6, we discuss the quantum mechanics of the oscillator 
system and the unitary transformations which correspond to canonical transformations in classical 
mechanics. In Sec. 7, we discuss physical applications of the formalism developed in this note. 
Finally, in Sec. 8, we discuss scale transformations in phase space and their implications in 
measurement theory. 


2 Linear Canonical Transformations in Classical Mechan- 
ics 


For a dynamical system consisting of two pairs of canonical variables x u pi and x 2 ,P 2 > we can 
introduce the four-dimensional coordinate system: 

071 , 172 . 73 . 94 ) = (zi,X 2 ,PuP2). ( 1 ) 

Then the transformation of the variables from rj x to is canonical if 


where 
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For linear canonical transformations, we can work with the group of four-by-four real matrices 
satisfying the condition of Eq.(2). This group is called the four-dimensional symplectic group 
or Sp(4). While there are many physical applications of this group, we are interested here in 
constructing the representations relevant to the study of two coupled harmonic oscillators. 

It is more convenient to discuss this group in terms of its generators G, defined as 


M = exp (— iaG), 


(3) 


where G represents a set of purely imaginary four-by-four matrices. The symplectic condition of 
Eq.(2) dictates that G be symmetric and anticommute with J or be antisymmetric and commute 
with J. 

In terms of the Pauli spin matrices and the two-by-two identity matrix, we can construct the 
following four antisymmetric matrices which commute with J of Eq.(2). 
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The following six symmetric generators anticommute with J. 
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These generators satisfy the commutation relations: 


(4) 


(5) 


\J * i J j] — ic** jk^kj [J , , K j] — icijkhky Kj] — [Qi, Qj\ itijkJk) 

[Jif Jo\ ~ 0 , [Ki,Qj] — iSijJ 0 , 

[Ji, Qj] = itijkQk, [Ki, Jo] = iQii [Qi, Jo] = -%Ki. (6) 


The group of homogeneous linear transformations with this closed set of generators is called the 
symplectic group Sp( 4). The J matrices are known to generate rotations while K and Q matrices 
generate squeezes [6] 

It is often more convenient to study the physics of four-dimensional phase space using the 
coordinate system 

= (zi,pi,x2,P2). (") 
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The transformation from (71,7^73,74) is 
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and the J matrix becomes 

/ 0 1 0 0\ 

, -1 0 0 0 

0.0 0 1 ■ 

0 0 - 10 / 

In this new coordinate system, the rotation generators take the form 


( 8 ) 


( 9 ) 



J ' = T 

f 0 cr 2 \ 

Ua 0 J ’ 
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/ 0 / ’ 
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(10) 

The squeeze generators become 
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(11) 


When we deal with canonical transformations of functions of the coordinate variables, we have 
to use the differential operators. The rotation generators are [6] 
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It was noted that there are two convenient coordinate systems in phase space, namely those of 
Eq.(l) and Eq.(7). The above differential forms are applicable to both coordinate systems. They 
of course satisfy the commutation relations given in Eq.(6). 

It is remarkable that these operators are also applicable to the Wigner phase-space distribution 
function which is constructed from the Schrodinger wave function [6]. It is also remarkable that 
there are unitary transformations on the wave function which lead to canonical transformations 
of the Wigner function in phase space [6]. 


3 Construction of the 0(3,2) deSitter Group by Solder- 
ing Two Lorentz Groups 

In Sec. 2, we constructed the ten generators of canonical transformations acting on two pairs of 
canonical variables. The mathematics is straight-forward, but it is not too comfortable to study 
physics with ten independent parameters. We can have a better physical picture if we can study 
the problem in terms of concrete physical examples with smaller symmetries. 

The deSitter group 0(3, 2) is known to be locally isomorphic to the group Sp( 4). Indeed, as we 
shall see in this section, the notations for the generators of Sp( 4) given in Sec. 2 are the natural 
notations for the deSitter group. Thus, one way to study 5p(4) is to study 0(3, 2). In this section, 
we shall study 0(3,2) by constructing it by soldering two 0(3,1) Lorentz groups. 

In the space-time of (x,y, z,t,s), where x,y,z are three space-like variables and t and s are 
two time-like variables, we can consider two 0(3, l)-like Lorentz groups in the spaces of (x, y, -O) 
and (x,y,z,s) respectively. The generators of rotations applicable to the three-dimensional space 
of x, y and z are 
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The Lorentz boosts in the subspace of (x, y , z, t ) are generated by 

/0 0 0 * 0\ /0 0 0 0 0\ /0 0 0 0 0\ 

00000 0 0 0 ■ i 0 0 0 0 0 0 

A'i = 0 0 0 0 0 , I<2 = 0 0 0 0 0 , I< 3 = 0 0 0 *• 0 . (15) 

* 0 0 0 0 Of 000 00 ? 00 

Vo 0 0 0 0 / Vo 0 0 0 0/ Vo 0 0 0 0 / 

These three boost generators, together with the rotation generators of Eq.(14), form a closed 
Lie algebra for the Lorentz group applicable to the four-dimensional Minkowski space of ( x , y, z, t). 
The same is true for the space of (x, y, z, s) with the boost generators: 

/ 0 0 0 0 * \ /0 0 0 0 0\ /OOOOOX 

00000 0000? [ 0 0 0 0 0 
Qi = 0 0 0 0 0, £ 2 = 0 0 0 0 0 , Q 3 = 0 0 0 0 ? . (16) 

00000 00000 00000 

V* 0 0 0 0/ Vo ? 0 0 0/ Vo 0 ?■ 0 0/ 

The above two Lorentz groups have nine generators. If we attempt to form a closed set of 
commutation relations, we end up with an additional generator 

/0 0 0 0 0 \ 

0 0 0 0 0 

Jo — 0 0 0 0 0 , (17) 

0 0 0 0 -? 

Vo 0 0 ?■ 0 / 

which will generate rotations in the two-dimensional space of s and t. These ten generators form 
a closed set of commutations relations. 

We started with two 0(3, 1) Lorentz groups. Each Lorentz group has its own rotation subgroup. 
In the present case, both Lorentz groups share the same rotation subgroup. This is how these two 
groups are soldered. 

It is remarkable that this set of commutation relations is identical to that of Eq.(6). The 
group 0(3,2) is locally isomorphic to Sp{ 4). The group 0(3,2) occupies a very important place 
in relativity and elementary particle physics simply because it contains two Lorentz groups as 
its subgroups. The local isomorphism between 0(3,2) and Sp(4) enables us to study this group 
in terms of linear canonical transformations in classical mechanics or in the Wigner phase-space 
picture of quantum mechanics. 

4 Construction of the Sp(4) Symmetry Group for Cou- 
pled Oscillators by Soldering two Sp(2) Groups 

For two uncoupled oscillators, we can start with the coordinate system: 

(^li&i&J^-l) = { x l,Pl,X2iP2) ■ (18) 
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Since the two oscillators are independent, it is possible to perform linear canonical transformations 
on each coordinate separately. The canonical transformation in the first coordinate system is 
generated by 


Ax 


1 / ( 7 2 0\ 

2 V 0 0 )' 


B U 


_£/< j 3 0 \ r i( 0 \ 

~ 2 V 0 0/ ’ 1 2 V 0 0 )' 


(19) 


These generators satisfy the well-known commutation relations: 


[j4x,i?i] = iCi, [B\,C\] — —Mi, 


[Ci, j4i] — iB\ . 


( 20 ) 


It is also well known that this set of commutation relations is identical to that for the (2 + 
l)-dimensional Lorentz group. Linear canonical transformations on the second coordinate are 
generated by 


A 2 




Bo = - 


(° °) 

Vo <r 3 / ’ 




( 21 ) 


These generators also satisfy the commutation relations of Eq.(20). We are interested here in con- 
structing the symmetry group for the coupled oscillators by soldering two Sp( 2) groups generated 
by A U B 1 , Ci and A 2 ,B 2 ,C 2 respectively. 

It will be more convenient to use the linear combinations: 


A+ — A\ + A 2 , 5+ — B\ + i?2, 

C + = Cl + c 2 , 



A_ = A x - A 2 , = fli - B 2j 

c_ = Cx - c 2 , 


(22) 

These matrices take the form 




*). :) 


V\ / 1 


*--;(? -J’ b -=H? -J 

I* C - = Ko 

0 )• 

(23) 

The sets (A+, B+,C+) and (A + , B-,C-) satisfy the commutation relations of Eq.(20). The 
is true for (A_,i? + ,C_) and (/!_,£_, C+). 

Next, let us couple the oscillators through a rotation generated by 

same 

*-s(! V) 



(24) 

Then, A 0 commutes with A+,5+,(7+, and the following commutation relations generate 
operators A 3 , B 3 and C 3 : 

new 

[Ao, A_] — iA 3 , [Aq, BJ\ — iB 3 , 

[Aq, £?_] — iB 3 , 


(25) 

where 

?)• ?) 


0 J 

( 26 ) 
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There are now ten generators. They form the closed set of commutation relations of Eq.(6), if we 
identify these matrices as 

A+ = - Jo, A_ = —Jz, A 3 = - Ji, A 0 — J 2 , 

B+ — K 2 , B- — —Q\, B 3 = Q3, 

= Q 2 , C- = I< x , C 3 = (27) 

where the J, and Q matrices are given in Eq.(10) and Eq.(ll). 

In this section, we started with the generators of the symmetry groups for two independent 
oscillators. They are A\, B X} C\ and Ai, B X) Ci. We then introduced another generator A 0 to solder 
them up. This processes produced three additional generators A3, £3, C 3 which are —J\, Q 3 ,and — 
K 3 respectively. It is remarkable that K 3 ,Q 3 and J 0 form the set of generators for another Sp( 2 ) 
group. They satisfy the commutation relations 

{Qz, K 3 \ = -i Jo, [A" 3 , Jo] = -iQz, [Qz, Jo] = iK 3 (28) 

This symmetry group will play the major role in decoupling the coupled oscillator problem. 


5 Reparametrization of Coupled Oscillators 


Let us consider a system of two coupled harmonic oscillators. The Hamiltonian for this system is 

H = — | — pjH p 2 A A! x\ + B*x 2 4- C / XiX2l . (29) 

2 l m 1 m 2 ) 

where 

A' >0, B l > 0, 4A'£' - C" 2 > 0. (30) 

By making scale changes of x x and x 2 to (mi/m 2 ) l B Xl and (m 2 /mi) l Bx 2 respectively, it is possible 
to make a canonical transformation of the above Hamiltonian to the form [14, 15] 


H = ~ (pj 4- p]} + ^ { Ax\ + Bx 2 + CxiXz} , 


(31) 


with m = (m 1 m 2 ) 1 ^. This transformation is generated We can decouple this Hamiltonian by 
making the coordinate transformation: 


fy\\_fcos(al2) - sin(a/2) W x 2 \ 
\ J/2 / V sin(a/2) cos(a/2) /\x 2 / 


Under this rotation, the kinetic energy portion of the Hamiltonian in Eq.(31) remains invariant. 
Thus we can achieve the decoupling by diagonalizing the potential energy. Indeed, the system 
becomes diagonal if the angle a becomes 


tan a = 


c 

B- A' 


(33) 


This diagonalization procedure is well known. 
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We now introduce the new parameters K and y defined as 


, A+B + J(A - B) 2 + C 2 

K = ,/AB-Cy 4, exp (-2 y) = ^ . 


(34) 


in addition to the rotation angle a. In terms of this new set of variables, A , B and C take the 
form 


A = K ( e 2 ^ cos 2 ^ + e 2r> sin 


B = K ^e 2T? sin 2 ^ + e 2,1 cos 2 

/I = K ( e“ 2r? — e 27? ) sin a. 
the Hamiltonian can be written as 

H = ^ {?i + ? 2 2 } + y ( e27, Pi + e 2tj ?/2 } > 

where y\ and yi are defined in Eq.(32), and 

qi \ _ ( cos(a/2) - sin(a/2) \ f pi 


*!)• 
!)■ 


92 


sin(a/2) cos(a/2) / \P 2 


W2)WPA 

a/2) ApJ' 


(35) 


(36) 


(37) 


This form will be our starting point. The above rotation together with that of Eq.(32) is generated 
by Jo- 

If we measure the coordinate variable in units of (ml\ )^ 4 , and use (m/i )“ li/4 for the momentum 
variables, the Hamiltonian takes the form 


H = A" (e % 2 + e v yl) + v (e n q\ + e ’’y 2 ) , 


(38) 


where u = ^ K jm. If y = 0, the system becomes decoupled, and the Hamiltonian becomes 


H = ^ (Pi + *i) + | (p 2 + x l) • ( 39 ) 

In Sec. 8, we will be dealing with the problem of what happens when no observations are made 
on the second coordinate. If the system is decoupled, as the above Hamiltonian indicates, the 
physics in the first coordinate is solely dictated by the Hamiltonian 

= | (p 2 + x i) • ( 4 °) 

It is important to note that the Hamiltonian of Eq.(39) cannot be obtained from Eq.(38) by 
canonical transformation. For this reason, the Hamiltonian of the form 


H f — 2 ( e 71 ^ + eV y^) + 2 ( eV ^ + e n yt) 

may play a useful role in our discussion. This Hamiltonian can be transformed into the decoupled 
form of Eq.(39) through a canonical transformation. 
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6 Quantum Mechanics of Coupled Oscillators 

It is remarkable that both the Hamiltonian H of Eq.(38) and H’ of Eq.(41) lead to the same 
Schrodinger wave function. If y\ and y 2 are measured in units of (mK) 1 ' 4 , the ground-state wave 
function for this oscillator system is 

0o(*i,2 2 ) = ^ exp {~\( eVy i + e_ ”y2)}- ( 42 ) 

The wave function is separable in the j/i and y 2 variables. However, for the variables Xi and x 2 , 
the story is quite different. If we write this wave function in terms of and x 2 , then 


*fi(xi,x 2 ) = -J=exp{-^ 


a , a 2 

e n (xx cos — — x 2 sin — ) 


.a or . ol 1 

+e v (xi sin - + x 2 cos - ) 1 1 . 


If r) = 0, this wave function becomes 


(43) 


^o(a?i,x 2 ) = -^exp|-^(xj + x 2 )|. (44) 

For other values of 77, the wave function of Eq.(43) can be obtained from the above expression by 
a unitary transformation. 

^ m2 (o, (x| )0 m2 (*^2)? (45) 

m 1 m2 

where ^ m {x) is the m Lh excited state wave function. The coefficients A mim2 (r}) satisfy the unitarity 
condition 


|4»im 2 ( a . , /)| a = L ( 46 ) 

m\ 7712 

It is possible to carry out a similar expansion in the case of excited states [16]. 

As for unitary transformations applicable to wave functions, let us go back the generators of 
canonical transformations in classical mechanics in Eq.(I2) and Eq.(13). As was stated before, they 
are also applicable to the Wigner phase-space distribution function. The canonical transformation 
of the Wigner function is translated into a unitary transformation of the Schrodinger wave function. 
There are therefore ten generators of unitary transformations applicable to Schrodinger wave 
functions [6, 4]. They are 


Ji 

J 3 

kx 


- ^(4 a 2 + 4 a i)> h = Yi (4 a 2 - 4 a i) > 

), Jo = ^ (4«i +^4) i 




— { a \ a l — °2 a 2 


^44 + a x a x - 44 —<*202) , 
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I<2 


I <3 

Qi 
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if ft It 

— ( <2j Gj — d\d\ “1" ^2^2 *” ^2^2 

^ (44 + a x a 2 J , 


)■ 


i 

4 

1 

4 


ajaj — a\Q\ — a[a\ + 02^2^ > 

a\a\ -f a x a x + a\a\ + a 2 «2^ ? 


Q 3 


2 [44 - °i a 2 


(47) 


where at and a are the step-up and step-down operators applicable to harmonic oscillator wave 
functions. The above operators also satisfy the commutation relations given in Eq.(6). 


7 Wigner Functions and Uncertainty Relations 

The Wigner phase-space picture of quantum mechanics is often more convenient for studying the 
uncertainty relations. Unitary transformations in the Schrodinger picture can be achieved through 
canonical transformations in phase space. It has been known that canonical transformations are 
uncertainty-preserving transformations. They are also entropy-preserving transformations [17]. 
Are there then non-canonical transformations in quantum mechanics? 

In his book on statistical mechanics [18], Feynman raises the issue of the rest of the universe 
in connection with the density matrix. Feynman divides the universe into two parts. We make 
measurements in the first part, but are not able to measure anything in the second part. The 
second part is Feynman’s rest of the universe. Indeed, the density matrix plays the essential role 
when we are not able to measure all the variables in quantum mechanics [19, 20]. 

In the present case of coupled harmonic oscillators, we assume that we are not able to measure 
the x 2 coordinate. It is often more convenient to use the Wigner phase-space distribution function 
to study the density matrix, especially when we want to study the uncertainty products in detail 

[!5, 18 ]. 

For two coordinate variables, the Wigner function is defined as [15] 

Vr(xj,x 2 ;pi,p 2 ) = (T J exp{-2z(p 1 j/ 1 +p 2 y 2 )} 

x t/>*(xi + 1/1 , x 2 + j/ 2 )V’(^i - Pi, *2 - V2)dyidy 2 . ( 48 ) 

The Wigner function corresponding to the oscillator wave function of Eq.(43) is 

/ 1 \ 2 f (% (x 

W(x u x 2 \pi,p 2 ) = - exp — e ,? (x 1 cos — — x 2 sin — ) 2 
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( 49 ) 


. a °\2 -ru ° ' Q \2 

— e v (x x sin — + x 2 cos — ) 2 — e (pi cos — — p 2 sin — ) 

. a o.o’l 

-e r '(pisin- + p 2 cos-) J . 

If we do not make observations in the x 2 p 2 coordinates, the Wigner function becomes 

W{x u pi) = J W{x u x 2 ,p u p 2 )dx 2 dp 2 . 

The evaluation of the integral leads to 


(50) 


2 7] sin 2 a) j 

_+ 4 

cosh 77 - sin rj cos a cosh jy 4- sin rj cos a 

This Wigner function gives an elliptic distribution in the phase space of xj and pi. This distribution 
gives the uncertainty product of 

(Ax) 2 (Ap) 2 = ^(1 + sinh 2 7/sin 2 a). (52) 

This expression becomes 1/4 if the oscillator system becomes uncoupled with cx 0. Because xj 
is coupled with x 2 , our ignorance about the x 2 coordinate, which in this case acts as Feynman’s 
rest of the universe, increases the uncertainty in the Xi world which, in Feynman’s words, is the 
system in which we are interested. 

In the Wigner phase— space picture, the uncertainty is measured in terms of the area in phase 
space where the Wigner function is sufficiently different from zero. According to the Wigner 
function for a thermally excited oscillator state, the temperature and entropy are also determined 
by the degree of the spread of the Wigner function phase space. 


x exp 



W(xi,x 2 ;pi,p 2 ) = 


x 2 ( 1 + sinh 


8 Scale Transformations in Phase Space 


In addition to the ten generators given in Eq.(10) and also in Eq.(ll), we can consider the scale 
transformation in which both the position and momentum of the first coordinate are expanded 
and those of the second coordinate contracted. The Hamiltonian given in Eq.(38) suggests such a 
transformation, and the transformation can be generated by 


So = 


if 1 

2 VO 



(53) 


This matrix generates scale transformations in phase space. The transformation leads to a radial 
expansion of the phase space of the first coordinate [21] and -contracts the phase space of the 
second coordinate. What is the physical significance of this operation? As we discussed in Sec. 
7, the expansion of phase space leads to an increase in uncertainty and entropy. Mathematically 
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speaking, the contraction of the second coordinate should cause a decrease in uncertainty and 
entropy. Can this happen? The answer is clearly No, because it will violate the uncertainty 
principle. This question will be addressed in future publications. 

In the meantime, let us study what happens when the matrix Sq is introduced into the set of 
matrices given in Eq.(10) and Eq.(ll). It commutes with Jqi *^ 3 ? ^ 1 ^ 2 ? Qi? an< ^ However, its 
commutators with the rest of the matrices produce four more generators: 

[So, -A] = 2 (aj 0 2 ) ’ (So, = 2 ( / o)' 

[*•*1-5(4 ?)• T)- < 54 > 

If we take into account the above five generators in addition to the sixteen generators of Sp( 4), 
there are fifteen generators. They form the closed set of commutation relations for the the group 
5’T(4,r). This 5L(4,r) symmetry of the coupled oscillator system may have interesting physical 
implications. 
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